Joint free energy distribution in the random directed polymer problem 
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We consider two configurations of a random directed polymer of length L confined to a plane 
and ending in two points separated by 2u. Defining the mean free energy F and the free energy 
difference F' of the two configurations, we determine the joint distribution function Vl,u{F , F') 
using the replica approach. We find that for large L and large negative free energies F, the joint 
distribution function factorizes into longitudinal {Pl.u{F)) and transverse {Pu{F')) components, 
which furthermore coincide with results obtained prevously via different independent routes. 

PACS numbers: 05.20.-y 75.10.Nr 74.25. Qt 61.41.+e 



Directed polymers subject to a random potential ex- 
hibit non-trivial behavior deriving from the interplay be- 
tween elasticity and disorder; numerous physical systems 
can be mapped onto this model and the topic has been 
the subject of intense studies Despite its undis- 

putable importance, our knowledge on this generic prob- 
lem is still limited. Traditionally, the main focus is 
on the free energy distribution function, for which two 
types of analytical solutions are known for the (1 + 1)- 
dimensional polymer confined to a plane (see Fig. 

[T]): one class addresses the 'longitudinal' problem and 
determines the distribution function Pl{F) for the free 
energy F of a polymer of length L and fixed endpoint 
y = [3, 0, 0, HI, 01 , while the other concentrates on the 
'transverse' problem aiming at the distribution function 
Vu{F') involving the free energy difference F' = F'^—F^ 
between two configurations with cndpoints at ?; = ±w 
0) B @]) assuming no dependence on the mean energy 
F ={F'^ + F~)/2 in the limit L^oo. Both approaches 
have been helpful in finding the wandering exponent C 
[lot of transverse fluctuations 5y{L) oc L'^ of the polymer. 
On the other hand, questions how the result for Vu{F') 
is approached from finite L and how the transverse and 
longitudinal problems are interrelated have remained un- 
clear; it is the purpose of this letter to shed light upon 
these issues. 

Here, we generalize the task of finding the free en- 
ergy distribution function for a polymer of length L by 
studying two configurations of the string ending in two 
points separated by 2u, see Fig. [1] and treating both 
the mean free energy F and the free energy difference 
F' as relevant variables. The two-point object F' relates 
to the natural variable appearing in the Burgers prob- 
lem [7j], while for u = the variable F reduces to the 
free energy F of a single configuration studied in Refs. 
d, H, 0, Hj @| • Our new scheme then should allow us 
to place the previous results for Vl{F) and Vu[F') into 
a common context. Using the replica approach, we de- 
termine the joint distribution function Vl^u^F, F') and 



prove (for a (5-correlated disorder potential) the separa- 
tion VlAF^F') = Pl,u{F) Pu{F') in the limit of large 
L and for large negative values of the mean free energy 
F . Furthermore, we derive the form of the two fac- 
tors Pl,u{F) and Pu{F'): on the one handjWe find that 
Pl,u{F) has the same form as Zhang's tail [3| for VLiF). 
On the other hand, to our surprise, we find that the trans- 
verse part Pu{F') exactly coincides with the stationary 
distribution function Vu{F') of the Burgers problem [3], 
although our solution is associated with rare events in 
the far-left tail, while the result of Ref. (3] describes an 
equilibrium situation reached in the limit L — > oo. In the 
following, we first describe the previous replica analysis 
leading to the distribution function Vl (F) and its poten- 
tial pitfalls and then proceed with the derivation of the 
joint distribution function Vl.u{F ,F'). 

We consider an elastic string (elasticity c) directed 
along the x-axis within an interval [—L, 0] and subject 
to a disorder potential driving the displacement 

field y{x), see Fig. [T] its energy is given by 

H[y{x)-V]= j"" ^dx[^[d^y{x)\^ + V[x,y{x)]]. (1) 




FIG. 1: Illustration of thermally averaged trajectories 
{y{x))th of a random directed polymer in a fixed disorder 
potential V{x,y). We let the polymer start in an arbitrary 
position at a; = ~L and fix the displacement y = m at a; = 0. 
Forcing the polymer to end in y — —u produces an alternative 
average trajectory on a distance ~ I. Our focus then is on the 
calculation of the joint distribution function Vl,u{F , F') for 
the mean free energy F — {F'^ + F~)/2 and the free energy 
difference F' = F^ — F~ of the two configurations. 
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The disorder average is carried out over a Gaussian distri- 
bution with zero mean {V{x,y))v = and a (5-correlator 
{V{x, y)Vix', y')). = Uq S{x - x')6iy - y'). 

The standard procedure ^ leading to the distribution 
function Vl (F) starts from the partition function (we set 
the Boltzmann constant ks equal to unity) 

/•a(o)=o 

Z{L-V)^ V[y{x)]eM~H[y{x)-V]/T), (2) 

Jy{-L)=0 

providing us with the free energy F{L; V) = —TlnZ. 
The A^-fold replication of the polymer and subsequent 
averaging over disorder realizations V maps the problem 
to N quantum bosons with local interactions —UoS{y — 
y')/T. In the large L limit, the replica partition func- 
tion Zr{N;L) = {Z^{L;V))v is dominated by the 
ground state of the quantum problem, which can be ob- 
tained from a Bethe Ansatz solution Zr{N;L) oc 
exp[-(A^(i^}^, + e3N^L)/T] (with eg = ~cU§/24T^). 
Exploiting the relation between the replica partition 
function Zr{N;L) and Vl{F) as given by the bilateral 
Laplace transform 



/oo 
dFVdF) expi-NF/T) 
-OO 



(3) 



allows one to show [3] that the far-left tail of the dis- 
tribution function VlIF) assumes the form VLiF) oc 
exp[— (2/3)(|F|/F»)'^/^] with the characteristic free en- 
ergy scale F^ = {cU^L/^T^f/^ cx L^s. 

In Ref. i, an attempt has been made to use the re- 
sult for Zr{N\L) and extract the third moment of the 
distribution function Vl{F). While predicting a wrong 
prefactor [13], this approach also misses to produce re- 
sults for other moments. The reason for this failure was 
identified by Medina and Kardar [l^, 14 1, who pointed 
out that the two limits L ^ oo (allowing to ignore excited 
states) and — > (providing the irreducible moments 
((F'=))v = {-T)''d'^\n{Zr{N; L)\N^a of the distribution 
function) do not commute. To obtain estimates for mo- 
ments, the assumption has to be made that the distribu- 
tion function Vl{F) is governed by a unique free energy 
scale cx L^^^; although this assumption cannot be 
expected to work for the very distant (non-equilibrium) 
tails, it turns out that its validity indeed extends to the 
far-left tail in the present problem, but does not for the 
far-right tail Summarizing, the original Bethe Ansatz 
solution allows one to find the (far-left) tail of the dis- 
tribution function [311 but cannot a priori provide infor- 
mation on its body [IsL [3] as this requires knowledge 
of the behavior of Zr{N; L) for TV ^ 0. 

Here, we study a different setup involving two con- 
figurations of a polymer with length L ending in points 
separated by 2u; we define the mean free energy F ~ 
(F+ -I- F-)/2 and difference F' = F+ - F- , with 
F^ = F{L, ±tt; V) the free energies of polymers ending 



in 2/(0) = ±M. The quantities F and F' are random vari- 
ables and we are aiming for the joint distribution func- 
tion Vl,u{F, F'). We define the replica partition function 
Z'^{n,m; L,u), which can be expressed as the bilateral 
Laplace transform of Vl,u{F,F'), 

Z[.{n, m; L, u) = {Z'\L, u; V)Z"^{L, -u; V))v (4) 

'iF+ rriF^ \ / {n+m)F (n~m)F' 
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-{n+m)F/T ^-(n-m)F' /2T 



The average over disorder realizations V provides us with 
the replica partition function in the form 



Z'^{n, m; L, u) 



- " rya(0)=u -1 r n+m „y^(o) = -u- 

n hiy^i^)] n / 

a=l •' J \-a=n+l •' 



xexp(-i/„+„[{2;a(a:)}]/r), (5) 
with the replica Hamiltonian 

/o " 

E S[ya{x)-y,{x)]}. (6) 

a, 6=1 

The replica partition function Z[.{n,m] L,u) describes a 
system with n+m trajectories ya{x) (a = 1, . . . , n+m) of 
which n traces terminate at the point u, while the other 
m trajectories end in the point —it; we adopt free initial 
conditions at x = —L as implied by the absence of 
any restriction on ya{—L) in ([S]). All these trajectories 
are coupled by the attractive potential —U^Slija — yb)/T 
deriving from the disorder correlator. 

We use the standard way ^ to map the path integral 
([5]) to a Schrodinger problem: allowing the n + m trajec- 
tories to end in an arbitrary point y = (?/i, . . . , ?/„+,„), 
we define the wave function ^'(y; x) — Z'^{n^ m; L + x, y) 
which satisfies the imaginary-time Schrodinger equation 
—Tdx'^{y",x) = H'^{y;x) with the initial condition 
\I'(y; — X) = 1. The Hamiltonian reads 



a, 6^1 



2 n+m 



Ur 



71 +m 



2c 



Vb) (7) 



and describes n+m particles of mass c/T^ interacting via 
the attractive two-body potential —UoS{y — y')/T. The 
partition function ([5]) is obtained by a particular choice 
of the final-point coordinates, Z^{n,m; L,u) — ^'(u;0) 
with u = (ui, . . . ,u„ = u;u„+i, . . . , u„+m = -"")• 

The expansion of ^'(y;a;) in terms of eigenfunctions 
^'if.a = cxp[iKY,aya/{n + m)]'0a ({i/a }) of ^ iuvolvcs 
a center of mass component and a factor ipaiiVa}) de- 
pending only on relative coordinates ya — yt- Our choice 
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of free initial condition ^'(y; —L) = 1 implies a vanish- 
ing center of mass momemtum A' = and our expansion 
assumes the simplified form 



(8) 



with Ea the eigenenergies. The coefficients Ca = 
(■0Q,|4'(— i) )/('!/'« I "00 ) follow from the initial condition 
^'(y; — L) = 1 with the scalar product {ip\<j>) = 
dya] SlEa Va/[n + m)] 
In the limit of large L, fixed u, and for integer n,m > 1 
(see below for a detailed discussion on limits and scaling u 
versus L) the sum in ([8]) is dominated by the ground state 
wave function ipo, for which the Bethe Ansatz provides 
the solution [11 1 



-00 (y) = exp(^-K^ |ya - yb\ 

a,b 



(9) 



with the inverse length k = cUq/AT^ and the energy [l3| 

Eo{n + m) = -cU^{n + m)[{n + - l]/2iT'^. (10) 

The normalization (^ol^o) — {n + m)/(4K)"+™^-^r(n + 
to) and the matrix element {ipol'i' L)) = {n + 
TO)/(2K)"+™-ir(n + to) provide the result ^'(y;0) = 
2„+m-ig-/3BoL^g(-y) j^jj^j evaluating (O at the endpoint 

u, we obtain the expression ipo{'^) = exp[— 4K|M|nTO] and 
hence 



Z'^{n, to; L, u) 



_ 2n+m — l ^ — Eoin+r, 



)L/T ^—iK\u\nm 



(11) 



Rewriting the exponent 4K|u|nTO = K|M|[(n + m)^ + (n — 
to)^], we can factorize Z'^{n, to; L, u) = Z^(n + m; L, u) 
Z~ {n — to; u) with 



Z- = 



K,\u\{n~m 



(12) 
(13) 



depending only on the variables n + m and n — m, cf. ([¥]). 
Hence we find that the 'transverse' problem described by 
Z~ {n — m;u) can be separated from the (mainly) 'lon- 
gitudinal' part encoded in Z^{n + m]L,u). This sepa- 
ration into transverse and longitudinal factors is a cen- 
tral element of our solution and tells us that the joint 
distribution function Vl.u{F', F) as defined in ^ factor- 
izes as wen, Vl,u{F,F')^ PlAF)Pu{F'). Corres_pond- 
ingly, we find that the distribution functions Pl,u{F) and 
Pu{F') are related to the factors Z+(n -I- m;L,u) and 
Z~{n — to; u) through the bilateral Laplace transforms 



Z^{n + to; L, u) 



Zj. (n ■ 



+ 00 

dFPL,u{F)e 

— OO 

/ + 00 
dF'P^{F')e' 
-OO 



-{n+m)F/T^ (^-^^^ 



We note that the above results could be derived for fixed 
initial conditions ya{—L) = as well, however, in this 



case the factorization appears only in the limit L — > oo. 
Also, the restriction to to, n > 1 limits the accessible 
values of F to large negative values and restricts the fac- 
torization of Vl,u{F, F') to the far-left tail in F. 

The expression ([T^ for Z~ has been derived for pos- 
itive integer n,m > 1 and large L; its dependence on 
n — m defines Z~ on all integers and simple inspection 
of (jlSp allows us to (uniquely) infer the final expression 
for the free energy distribution function 



P 



/ r \i/2 / TF'^ \ 



Formally, the result can be obtained via analytic 
continuation of Z^ into the complex plane and use of the 
inverse Laplace transform (we define ^_ = (n — m)/2T) 



Pu[F')= ^Z-{2T^^;u)expi^^F'), 



(17) 



requiring an analytic continuation of Z^ to the imagi- 
nary axis. This procedure leads to the identical result 
(fTHll . however, without solid control on the analytic con- 
tinuation. The result (fTO)) coincides with the Gaussian 
distribution function for the velocities in the correspond- 
ing Burgers problem Q, including all numericals. This 
comes as a surprise and may indicate that the factoriza- 
tion, which we can prove for the far-left tail, may actually 
prevail throughout all values of F. 

Next, we analyze what information on Pl,u{F) can be 
extracted from Z^in + to; L, u). For u = 0, the distribu- 
tion function Pl,u{F) coincides with Vl{F), F = F, and 
the partition function Z^{n + m]L,u) with the ground 
state approximation of Zr{N; L), N — n + m. The par- 
tition function Z^ as given by (|12p is valid for positive 
N = n + m and provides, via (jl4p . information on large 
negative free energies F = F, i.e., the left tail of the dis- 
tribution function, Pl{F) cx exp[-(2/3)(|F|/F*)3/2] as 
calculated by Zhang 3]. Inserting this result back into 
Eq. and evaluating the integral via the method of 
steepest descent, one finds that the main contribution to 
the integral arises from values F ^ —{F^/T'^){n + m)'^; 
negative free energies such that —F > F^/T'^ then 
correspond to positive values n + m for which we can 
trust the expression for Z^(n + to; L, 0) and hence for 
Vl{F\ Going to finite u, we still can trust our result 
for Z+(n -I- m;L,u) provided that Fei_= cv?/2L < |F| 
(see below) and we find a factor Pl,u{F) of basically the 
same form as for u = 0. 

In order to assess the regime of validity of our re- 
sults, we have to study the contribution to Eq. ([5]) of 
excited states. For u = 0, the relevant excited state is 
the one with lowest energy; this state is one-fold ionized 
[lit ^^'^ its excitation 'energy' is given by AEL/T — 
(i^*/T)3(n + TO)(n + TO-l). With {n + m)^ ^_\F\T^ / F^ , 
ground state dominance then requires that F ^ T and 
combining this condition with the one above we find that 
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\F\ ^ max[F^/T'^,T]. Introducing the temperature de- 
pendent Larkin length L^iT) ~ T^/cU^, see Ref. El, 
this condition assumes the form |F| 3> max[T, TL/LJ- 
For large u, the most dangerous excited state involves 
two free clusters with n and m bound particles; with 
an 'energy' EL/T = -{F^/Tf[n(n'^ - 1) + - 
l)]/3 and no tunneling suppression through the excited 
state wave function, we find a difference in exponents 
[(i^*/T)'^(n + m) — 4:K\u\\nm, from which we obtain the 
condition \F\ ^ u^c/L] the combination with the re- 
strictions obtained before produces the overall condition 
\F\ > max[cu2/L,r, TL/LJ. Hence, for L > ic, typi- 
cal excursions Sy{L) oc L^/'^ are well within the domain 
of applicability of our results. 

In analogy with (jl7p . one might directly apply the in- 
verse Laplace transform to the approximate result 
as given by (|12p . Dropping terms linear in n -f m 
and choosing u — 0, one easily recognizes the integral 
representation of the Airy function, C~^[Z^{n + m)] oc 
Ai{-F/F^). The asymptotics at negative F — —\F\ of 
the Airy function agrees with Zhang's tail of the distribu- 
tion function Vl{F), as already noted above. However, 
pushing the free energy F to positive values, the charac- 
teristic oscillations of the Airy function are incompatible 
with the positivity of the distribution function Vl{F). 

Although the above simplified approach correctly ac- 
counts for the center of mass (COM) degrees of freedom, 
it still fails to produce a consistent result for Vt,{F). 
This observation is in line with a previous study [20|, 
where the COM motion was accounted for and a nega- 
tive mean square displacement {{Sy'^(L)) th)v was found 
in the N ^ limit, but contradicts to the claim in Ref. 
[lit that the inclusion of the COM motion leads to a 
consistent result. We attribute the severe problems ap- 
pearing in the derivation of VhiF) to the impossibility 
to analytically continue the ground state approximation 
of Zr{N;L) derived for integer > 1 and large L to 
values A' < 1: at A^ = 1, all the spectrum describing 
the relative motion between bosons collapses to and 
the former ground state energy reappears at A^ < 1 with 
positive energy, cf. PTI)) . As a result, there is no control 
on the relevant excitations in the regime A < 1. 

While the inconsistencies in the analytical continuation 
of the replica number N = n + m across unity are quite 
prominent in the longitudinal problem of finding Vl (F) , 
they appear much more subtle in the analogous calcula- 
tion of the transverse distribution function Vu{F'): Fol- 
lowing Ref. and setting n + m — in (j4]), the inte- 
gration over F could be trivially done and the inverse 
Laplace transform of — Z~ /2, see Eq. p3|) . provides 
a result for Vu{F') which, surprisingly, is correct up to 
a prefactor 1/2. Although the missing excitations entail 
merely a spoiled normalization in this case, the conse- 
quences of dropping excitations are much more drastic 
when dealing with fixed initial conditions where the pre- 
factor diverges. Since our above analysis oi Vl.u{F, F') 



does not rely on the replacement of n -I- m by zero it is 
devoid of these problems. 

In conclusion, we have calculated the joint distribution 
function Vh.uiF , F') for two polymer configurations with 
cndpoints separated by 2u, allowing us to discuss the 
longitudinal and transverse problems on an equal foot- 
ing, cf. (III). Starting from a modified replica approach, 
we make use of the Bethe Ansatz solution of the asso- 
ciated quantum boson problem: We find separability of 
the longitudinal and transverse problems at large lengths 
i, a transverse factor which, to our surprise, coincides 
with the stationary distribution in the Burgers problem 
0, and a longitudinal factor which agrees with Zhang's 
tail [s*!. The validity of these results is limited to large 
negative values of F ^ a consequence of keeping only the 
ground state wave function in the solution of the quan- 
tum problem. For a finite-width random potential corre- 
lator these conclusions remain (approximately) valid at 
not too low temperatures and not too large —F, whereas 
the decrease in temperature or the increase in —F lead 
to the disappearance of the factorization in 'Pl,u{F, F'). 
Further progress, particularly with respect to the longi- 
tudinal problem, seems to rely on a better understanding 
of the spectral properties of the quantum-boson problem 
in the regime < n + m < 1. 
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